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In cellular automata models a glider gun is an oscillating pattern of non-quiescent states that 
periodically emits traveling localizations (gliders). The glider streams can be combined to construct 
functionally complete systems of logical gates and thus realize universal computation. The glider gun 
is the only means of ensuring the negation operation without additional external input and therefore 
is an essential component of a collision-based computing circuit. We demonstrate the existence of 
glider gun like structures in both experimental and numerical studies of an excitable chemical system 
- the light-sensitive Belousov-Zhabotinsky reaction. These discoveries could provide the basis for 
future designs of collision-based reaction-diffusion computers. 
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I. INTRODUCTION 



The field of unconventional computing involves the design of novel computing devices based on principles of infor- 
mation processing in chemical, physical and biological systems [IlEllS]. Experimental prototypes include reaction- 
diffusion processors [H O O |71 [HI IHJ EHl IHl US! US] , extended analog computers [TJ , micro-fluidic circuits [IS] , and 
Plasmodium computers fl6! [T7] . These computing substrates utilize propagating patterns that encode for informa- 
tion, e.g., chemical waves in reaction-diffusion processors. Every micro- volume can act as both logical gate and wire 
meaning computation can be performed simultaneously at many points in the reactor. Therefore, it is counterintuitive 
to impose significant levels of geometrical constraint such as channels to direct information transfer. Instead, the aim 
is to utilize the natural dynamics of the system to implement computation. The theory of collision-based computing 
provides a rigorous computational framework to enable this type of implementation |18j . 

Collision-based computing was developed from pioneering results on the computational universality of the Game 
of Life [in], the billiard-ball model and conservative logic [20] and their implementation in cellular-automaton (CA) 
models [2T]. 

Collision-based computers employ compact mobile patterns (localizations) to encode elementary units of informa- 
tion. Truth values of logical variables are given by the presence or absence of localizations or alteration of localization 
type. The localizations travel in space and implement logical gates when they collide. There is no imposed architec- 
ture such as stationary wires, but instead trajectories of localizations constitute transient wires. The benefits are a 
hugely adaptable, massively parallel computing medium. 

A wide range of physical media are capable of supporting mobile localizations. These include solitons, defects 
in tubulin microtubules, excitons in Scheibe aggregates and breathers in polymer chains [18j . Theoretically, these 
phenomena should be candidates for implementing collision-based circuits, but in practice they are not easily observed 
or controlled. A relatively amenable experimental substrate based on the light-sensitive Belousov-Zhabotinsky (BZ) 
medium [35] was shown to be capable of implementing collision-based gates [531 Ull HSl [2S]- The sub-excitable 
phase of this reaction supports short-lived propagating localizations in the form of unbounded wave fragments that 
behave like quasi-particles. The reaction can be made sub-excitable by careful control of the light level. In the light- 
sensitive BZ reaction increased light decreases excitability by producing Br- ions, the main inhibitor of the reaction. 
The sub-excitable limit is between the excitable (plane waves propagate) and non-excitable (reaction insensitive to 
perturbations) phases of the reaction. 
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In previous experimental and numerical studies |23l |231 US [55] the design of functionally complete logical gates 
in BZ reactions was demonstrated. An outstanding problem from these studies was how to implement a negation 
operator without artificially feeding additional wave fragments into the computing medium. In CA models of collision- 
based computers (see [181), negation is implemented using autonomous generators of gliders, known as glider guns. 
This is done by colliding a stream of gliders with other traveling localizations representing data. Until recently there 
was no experimental proof that glider gun analogs existed in BZ systems. This paper presents new results showing 
the existence of glider gun analogs in both experimental and numerical studies of the BZ reaction. We also show that 
the proposed mechanism of glider gun formation in certain CA systems has close parallels with the formation of glider 
guns in real chemical systems. 

The paper is structured as follows. In Sect. |ll]we summarize results of glider gun formation in two distinct CA 
models and propose a mechanism of glider generation by glider guns. In Sect. |III| we describe the experimental setup 
and numerical parameters used to implement glider guns in the BZ reaction. In Sect. |IV| we discuss the process of 
glider gun formation in heterogeneous BZ systems. Sect. |V] describes several examples of glider guns implemented in 
numerical studies of the BZ reaction, including a computational interpretation of the results. The paper finishes with 
a section on approaches to obtaining glider guns in architectureless homogeneous BZ systems (Sect. |VI[ ). 

II. GLIDER GUN FORMATION IN CELLULAR AUTOMATA MODELS 

We propose that gliders in the following CA models are produced from interactions of spiral wave-like patterns 
that constitute the body of the gun and that these are analogous to wave fragments in a chemical reaction-diffusion 
medium. A comprehensive proof of this proposition is still somewhat distant as it would require exhaustive analysis of 
glider guns formed by all possible cell-state transition rules. The CA lattices 2+-medium [27 and Spiral Rule [551 US] 
automata demonstrate our reasons for suggesting this mechanism. 

A. Cellular automaton model of sub-excitable chemical media 

The 2+-medium cellular automaton is a two-dimensional orthogonal array of cells where each cell has three possible 
states: resting, excited and refractory. A cell's state is updated in discrete time depending on the states of its closest 
neighbors. A resting cell becomes excited if two neighboring cells are in the excited state. Excited cells become 
refractory and refractory cells adopt a resting state at the next time step . The 2"'"-medium is an excitable discrete 
system capable of universal computation j27j . In previous work, we implemented few-bit adders and multipliers in 
the 2+-medium [30l[31]. 

The simplest mobile localizations (gliders) in the 2+-medium are comprised of four non-resting states, a head with 
two excited states and a tail with two refractory states (2+-particles). Several mobile glider guns (where the body of 
the gun moves) that emit 2+-particles have been identified. To date, no stationary glider guns have been discovered 
in the 2"'"-medium. A mobile glider gun moving Westwards emitting 2"''-particles in the opposite direction is shown 
in Fig.jlja). 

Whilst traveling, the gun undergoes structural transformations (Fig.jljb)). On closer inspection the transformations 
involve the interaction of two spiral excitation waves within the body of the gun. These spirals collide to form the 
2+-particle (Fig. [ijb-c), direction of rotation shown by arrows). 

B. Reaction-diffusion based CA model 

The Spiral Rule CA |2H1[2S], derived from Wuensche's beehive rule [32], consists of a hexagonal array of cells, where 
each cell takes one of three states {2, 1,0}. Cells update their state depending on the state of their neighbors. The 
cell-state transition rules j28j can be represented as a matrix where i is the number of neighbors in state 2 and j is 
the number of neighbors in state 1. Therefore, according to its neighborhood a cell will update to cell-state Af[i][j] 
at the next time step (Fig. |2]). 

In [5S], the Spiral Rule CA was interpreted in terms of an abstract reaction-diffusion medium where '1' is the 
activator A, '2' is the inhibitor I, and '0' is the substrate S. Therefore, the transition matrix (Fig. [2| represents an 
abstract chemical reaction. Full details of the quasi-chemical reactions are described in [29] . 

As shown in Fig. [sja), a remarkable feature of the Spiral Rule CA is that it exhibits rotating glider guns. The gun 
is a spiral wave fragment with a few activator states, mainly concentrated at the head of the wave. We hypothesize 
that the tail of the spiral wave, containing one activator state, periodically loses stability and becomes detached. This 



lb) (c) 

FIG. 1: Glider gun in the 2"'"-medium. (a) The gun travels West and emits gliders (self-localized excitations) traveling East, 
(b) Oscillation cycle of the gun. (c) Mechanics of glider generation, p is the gun's velocity vector, g is the direction of the glider 
stream, r shows the inward rotation of the wave tips. A solid circle indicates the excited state, an open circle indicates the 
refractory state and a dot indicates the resting state. 

process forms a propagating localization (glider) with one activator state and three inhibitor states (Fig. [sj^b)). The 
tail is regenerated and the mechanism repeated, resulting in streams of gliders. 

III. GLIDER GUNS IN THE LIGHT-SENSITIVE BZ REACTION 

A. Experimental 

A light-sensitive BZ reaction [35^ was used, as detailed in Fig. |4j A controllable network of excitable (low light 
0.394 mW-cm^^) and non-cxcitablc (high light 9.97 mW-cm~^) cells arranged in a checkerboard pattern was projected 
onto the gel surface ^35, ,36» i32j . 

B. Numerical 

The Oregonator model modified to account for the photochemistry was used |38j . 
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Model parameters u and v are the dimensionless concentrations of the activator and the oxidized form of the 
catalyst, respectively, q = 0.0002 and e are scaling parameters, / is a stoichiometric coefficient. = 1 is the 
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FIG. 2: Transition matrix for the Spiral Rule CA [25] , 

diffusion coefficient of u, while Dy = since v is immobilized in the gel. The system was integrated using the Euler 
method with a five-node Laplacian operator (V^), time step At = 0.001 and grid point spacing Ax = 0.62. is the 
rate of production of the photo-induced inhibitor. Low excitability corresponds to high values of <P (equivalent to 
high light) and vice versa. At = the medium is oscillatory. The medium is excitable up to <P = 0.034, weakly- 
excitable between <l> — 0.034 and <P — 0.0355 (the sub-excitable threshold) and non-excitable above — 0.0355. The 
model network consisted of alternating excitable {'I'l) and non-excitable ((f^) squares. Waves were initiated from an 
oscillating square (<? = 0) and directed into excitable squares of the network. 

IV. MECHANISM OF GLIDER GUN FORMATION IN BZ SYSTEMS 

Glider guns in BZ systems are rotating spiral waves which become trapped within single excitable cells in a network 
of light (non-excitable) and dark (excitable) cells. Excitable cells have up to four neighboring excitable cells joined 
by corner to corner junctions. Therefore, streams of gliders (excitation waves) can be emitted in four directions. 
Examples of glider guns in numerical and experimental BZ systems that emit glider streams in all possible directions 
are shown in Fig. [5j 

Glider guns are formed when spiral waves with trajectories pinned around non-excitable cells degenerate, giving 
spirals trapped in excitable cells. If waves are input to networks with high excitability fragmented stable waves are 
formed that travel through all corner to corner junctions of the excitable cells. However, if network excitability is 
reduced (#j and moved closer to the sub-excitable limit) spirals form as the fragmented waves become disconnected 
due to junction failure. These spirals are often complex interacting spirals with their trajectories pinned around non- 
excitable cells or regions. Therefore, fragments are incident repeatedly at certain junctions and eventually only critical 
amounts of activator (on the boundary between expansion and annihilation) can cross these junctions. These critical 
excitation waves form as the refractory period of the junction is only fractionally longer than the period of rotation 
of the original spiral(s). Critical excitation waves do not expand rapidly but instead drift across the cell. When they 
eventually expand they possess high curvature enabling them to curl back to form spirals trapped within excitable 
cells. Wave curvature increases with excitability [35] and the core of a spiral wave decreases |1D]. Therefore, excitable 
cells in the network must possess a critical size and excitability in order for these spirals to form. If excitability is too 
high the initial waves are stable and no spirals form in the network. If excitability is too low, spirals cannot form in 
excitable cells because wave curvature is too low. 

Spirals trapped within excitable cells have a shorter period of rotation than the parent spirals. Over time a number 
of these spirals form at various points in the network. This is because these spirals emit a proliferation of wave 
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FIG. 3: Glider gun in the Spiral Rule hexagonal CA: (a)-(f) oscillation cycle of the gun, (g) mechanics of glider production, 
g is the velocity vector of a glider traveling North-East, r indicates the rotation of the spiral wave, e is an excitation force. 
Activator state is shown by a solid circle, inhibitor state by an open circle and substrate by a dot. 



fragments meaning that the original mechanism of formation has an increased chance of being repeated. 

The glider guns formed in BZ systems exhibit strong parallels with glider guns formed in the Spiral Rule CA. In this 
CA a rotating spiral wave emits glider streams in all possible directions (Fig. [sj^a)). In the BZ system the interaction 
between the rotating spiral wave, the cell boundaries and junctions controls the firing of the gun. In Fig. |5] the 
junctions are relatively excitable and therefore, when a glider gun forms, gliders are emitted in all possible directions. 
So despite the periodic motion of the spiral, a super-critical amount of activator always reaches the junction and 
stable gliders form. 



6 




J 

FIG. 4: Light-sensitive catalyst (Ru(bpy)3^) [0.004M] loaded silica gel (D) ( 34 ) immersed in thermostated (G) Petri dish 
(E) containing catalyst-free BZ reagents [NaBrO3]=0.36 M, [CH2(COOH)2]=0.0825 M, [H2SO4l=0.18 M, [BrMA]=0.165 M. 
At this solution composition the reaction is oscillatory in the dark. A peristaltic pump (F) continuously feeds the reactor with 
solution and removes effluent (J). The reaction solution (H) is in an ice bath. The checkerboard is generated by computer 
(A) and projected (B) through a 455 nm narrow bandpass interference filter, lens pair and mirror assembly (C). Images are 
captured by a digital camera (I) and processed to identify chemical wave activity. 




(a) (b) 

FIG. 5: Glider guns formed in (a) experimental and (b) numerical studies of the BZ reaction. Grey arrows show the propagation 
direction of glider streams. All the streams of gliders come from a single glider gun (spiral rotating within one excitable cell) 
marked by a circle. Conditions (a) $1 = 0.394 mW-cm'^ and $2 = 9.97 mW-cm'^ (b) <Pi = 0.003, -Sg = 0.0255, e = 0.25, 
q = 0.0002 / — 1.11, Aa; = 1.3. Grid size (a) approximately 22x22mm, (b) 135x135 simulation points. 



In experiment, if excitability levels are reduced further then wave transfer is generally unstable because of added 
heterogeneity (gel imperfections, light dispersion etc.) when compared to numerical studies. In numerical studies 
we were able to map all excitability levels where glider guns formed. Glider guns that only emit wave fragments in 
selected directions were observed at lower excitability levels. In some cases glider guns fire intermittently into certain 
neighboring cells. This intermittent firing pattern can be explained by the periodic shift in rotation of the spiral 
(Fig.[6|. These 'unstable' glider guns are easier to fit into collision-based computational schemes as stable glider guns 
completely dominate network dynamics (see Sect. |v]). However, the fact that both types exist in the model indicates 
that unstable types may be observed in experiments where heterogeneity can be reduced. 
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FIG. 6: Activator concentration measured in numerical studies (at two opposite junctions in a cell) showing the formation of a 
glider gun. In region (I) fragments incident at the junctions are widely spaced as spiral waves in the network have trajectories 
pinned around non-excitable cells. In region (II) a spiral rotating in an excitable cell (glider gun) forms and closely spaced 
fragments are incident at the junctions. The amount of activator at each corner varies periodically because the spiral tip 
meanders within the cell. This observed periodicity controls the number and direction of glider streams emitted from the cell. 
Figures show dimensionless time and concentration units. 



V. COMPUTING WITH GLIDER GUNS IN BZ MEDIA 

In our numerical studies we have observed glider guns with 2, 3 and 4 stable streams of gliders (Fig.|7][8|. In addition, 
many glider guns exhibit an unstable intermittent glider stream from one or more of the 'empty' junctions (Fig.[7]|9|. 
The factors controlling the number of stable and intermittent glider streams are differences in the excitability of the 
cells and the complex periodic motion of the spirals trapped within excitable cells (Fig. [6]). 

As mentioned, in CA models gliders play an important role in implementing negation gates and thus ensuring 
computational universality. Glider guns formed in the BZ reaction can be used for the same purpose. The following 
reaction schemes show how 'data' wave fragments interact with glider streams in numerical studies of the BZ reaction. 
The work also shows how glider guns and glider streams can be manipulated and annihilated and the results interpreted 
in terms of computation. 



Fig. 10 shows the space-time evolution dynamics of two glider guns marked ei and 62 and auxiliary spiral wave 



fragments marked oi and 02 in Fig. 10 d). The spirals ei and 62 generate localized wave fragments (analogs of 
gliders) traveling North-East (stream gi) and North- West (stream 52)- At this network excitability gliders can only 
cross junctions that are opposite their direction of travel. Compare with Fig. [5] where excitability is higher and 



fragments emanating from the original glider gun cross all junctions in the network. In Fig. 10 the auxiliary spiral 
waves ai, a2 collide with the glider gun streams and cannot expand further into the network. The glider guns generate 
a constant output regardless of the presence of auxiliary 'data' wave fragments. Interestingly, where glider streams 
intersect they annihilate. The glider streams segregate the network into zones which is useful from a computational 
standpoint. This is because multiple reactions with 'data' fragments can be realized in parallel. 

In CA models of glider guns [29], the number of glider streams can be controlled by placing stationary localizations 
(termed 'eaters') at various points in the network. These 'eaters' destroy gliders by colliding with them. In the BZ 
medium we can employ localized spiral waves to annihilate redundant streams, as in Fig. [7] Here, the glider gun e in 
Fig. ^d) rotates counterclockwise generating three stable glider streams. The fourth stream (54) is annihilated by a 
spiral wave (a) rotating around a non-excitable cell. 

When constructing dynamically reconfigurable circuits it is important to be able to remove obsolete glider guns 
and generate or reposition existing glider guns by using other mobile localizations. This is useful where the bit-range 



(a) (b) (c) 




(e) :(f) (g) 

FIG. 7: (a)-(c) Excitation dynamics from numerical simulation of the BZ reaction showing cancelation of glider stream, 
(d) Reaction scheme showing that glider gun e generates four glider streams g\ - ■ ■ g^. The unstable stream is canceled via 
collision with auxiliary spiral wave a. (e) Following on from (c) the excitation dynamics show the formation of a second glider 
gun 62 (marked by a circle), (f) Cancelation of the parent glider gun ei by the daughter glider gun 62 (Reaction scheme shown 
in (g)). Network excitability (^1 = 0.0245 and $2 = 0.0465. Model parameters e = 0.11, Ax = 0.62. Grid size 180x180 
simulation points. 

of a circuit should be increased or when switching an arithmetic circuit between an adder and a multipher. This type 
of operation has been successfully implemented in CA models [251 [551 [31] . In CA models, glider guns are produced 
in collisions between several gliders. Glider guns can be extinguished via the collision with appropriate timing of a 
glider at a certain angle. In numerical studies of the BZ reaction we observed the same behavior. Initially, auxiliary 
spiral fragments and glider streams from ei combine to form a glider gun 62 (marked by a circle in Fig. |7]^e)). The 
'daughter' glider gun 62 is positioned in a such a way that stream (721 eventually annihilates the parent glider gun ei 
(Fig.jTjf-g)). 

Glider streams can be switched on and off periodically (blocked and unblocked) if two multi-streamed glider guns 
adopt a specific spatial arrangement (Fig. [8|. Each of the glider guns (61,62) would emit four streams of wave 
fragments but stream 512 stops all gliders traveling in a South-East direction from 62 (Fig. [8ja)(f)). Glider stream 
gii is periodically interrupted via interaction with glider stream 321- Glider gun 62 is continually firing but the 
glider stream 321 is temporarily interrupted at a point just below the gun due to interaction with glider stream 1712 
(Fig. |8]^c)). This gap in stream g^i allows stream 511 to be switched back ON (Fig. [8]^d),(e)). However, glider stream 
gii is switched off again for 360 dimensionless time units (t.u.) due to interaction with glider stream (721- This 
switching sequence was observed over many complete cycles. This type of switching circuit observed in the chemical 
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FIG. 8: (a)-(e) Excitation dynamics in numerical simulations of the BZ reaction showing the switching cycle between two 
interacting glider streams, (f) Schematic representation showing that glider gun ei initially generates four streams of gliders 
gii ■ ■ ■ g\4, (see (a)). Glider stream gw is temporarily switched off via interaction with stream (721 from glider gun 62 (see(b)). 
Glider stream gn is switched back on again after stream g2\ is interrupted by interaction with fragments from stream g\2 (see 
(c),(d)). Network excitability e^i = 0.0245 and $2 = 0.0465. Model parameters e = 0.11, Aa; = 0.62. Grid size 180x180 
simulation points. 



system provides a good example of how complex computational architectures may be achieved. 

The switching sequence of wave fragments in streams gn and (721 at the point of collision is described as follows: 



gii = 0_^ 0_^ . . . 

180 t.u. 360 t.u. 180 t.u. 360 t.u. 
321 = 01... 11... 101... 11... 1 

Fig. [9] shows how simple addressable/rewritable memory units can be achieved in numerical studies of the BZ 
reaction. In this scheme the reactor becomes segregated into two zones by the formation of glider gun e. In the upper 
zone a single auxiliary spiral a maintains a stable trajectory (Fig. [9][a)). This spiral is a simple analog of memory as 
it represents past interactions in the reactor. When fragments from the glider gun enter the cell marked with a circle 
(Fig.jojb)), spiral a is annihilated and the memory is erased (Fig. [oj^c)). This occurs periodically due to the complex 
rotation of the spiral (glider gun) causing changes in activator concentration at adjacent junctions. Therefore, once 
the original spiral has been erased this 'memory' zone of the reactor is periodically updated. Eventually, another glider 
gun forms (Fig. |9]^d)) which provides a constant input to the 'memory'. This reaction scheme shows how unstable 
junction transition causes a proliferation of glider guns. 



VI. GLIDER GUN-LIKE STRUCTURES IN HOMOGENEOUS BZ SYSTEMS 

The work described in the preceding sections involved the formation of glider gun-like structures in heterogeneous 
BZ networks with high and low excitability regions. Although these structures represent a step forward in terms of 
implementing a complete set of collision-based gates a heterogeneous architecture is imposed. 
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FIG. 9: Implementing simple memory analogs in numerical studies of the BZ reaction, (a)-(d) Excitation dynamics showing 
that stable spiral wave a (a) is annihilated by fragments emanating from the intermittent firing of glider gun e (b). This 
can be interpreted as the stored information being erased thus emptying the memory (c). Another glider gun is formed (d) 
representing a constant input to the memory (reducing the size of the addressable zone) . (e) shows the original reaction scheme 
where a represents the initial information stored in the memory. The memory unit is formed by the segregation of the reactor 
by the glider gun e and stream of gliders g. Network excitability $i = 0.0245 and — 0.0465. Model parameters e — 0.11, 
Ax = 0.62. Grid size 180x180 simulation points. 



Our previous studies of collision-based computing in experimental and numerical studies have utilized homogeneous 
BZ systems [231 Ull 113 121] • A homogeneous system has uniform excitability. In the case of the light-sensitive BZ 
reaction this is due to the projection of a uniform light intensity. If this uniform light intensity is close to the 
sub-excitable limit = 0.0355 model, « 5 mW-cm^^ experiment) then these homogeneous BZ systems exhibit 
traveling unbounded wave fragments. However, constant generators of these fragments have not been identified. Wave 
fragments are simply formed directly from waves injected into the homogeneous sub-excitable media. These wave 
fragments maintain stability over relatively long distances as at this excitability level their free ends do not curl to 
form stable spirals. 

One possible route to obtaining streams of mobile localizations in homogeneous BZ systems involves injecting 
fragments into a region of uniform light intensity set at a point just above the sub-excitable limit. In experimental 
studies of this type, we observed the formation of an unstable spiral that underwent spontaneous and periodic break- 
up. As the spiral rotated and expanded, the distance from the stable core increased and eventually the expanding arm 
of the spiral wave broke in two places producing wave fragments moving in a North-East direction (Fig. 11 a)-(c)). 
The continuous rotation and breakup of the spiral produced a localized stream of wave fragments. This experimental 
result exhibits strong parallels with the CA model described in Fig. [3] where the tail of an unstable rotating spiral 
spontaneously breaks and reforms to produce localized streams of fragments. However, this experiment involves 
injecting wave fragments into a weakly-excitable homogeneous medium and reliance on a highly non-linear process in 
the formation of an unstable spiral wave. 

To successfully implement collision-based computing schemes it is necessary to have a high degree of control over 
the spatial evolution of excitation wave fragments. Fig. |ll[ d)-(f) shows the same experimental setup with the addition 
of a dark region projected onto the otherwise homogeneous BZ system. Spontaneous initiation of target waves occurs 
in the dark region of the reactor. However, as waves move into the weakly-excitable region they spontaneously break 
to form a continuous stream of mobile wave fragments. This scheme imposes minimal architecture in the form of the 
dark area and offers benefits in terms of designing complex computing schemes as generators can easily be switched 
on and off. The position and firing frequencies of one or multiple generators can also be controlled. 





FIG. 10: (a)-(c) Excitation dynamics in numerical simulations of the BZ medium. Network excitability #i = 0.0245 and 
$2 = 0.0465. e = 0.11, Aa; = 0.62 (d) Schematic of glider gun interactions: ei and 62 are spiral waves (glider guns) generating 
glider streams gi and 32 that interact with auxiliary spiral waves ai and 02- Grid size 198x198 simulation points. 




FIG. 11: Gliders in a homogenous sub-cxcitable medium, (a)-(c) snapshots of excitation dynamics with spiral initiation, (d)-(f) 
snapshots of excitation dynamics with initiation via an oscillating dark region {$ = 0.034 mW-cm~^). Sub-excitable light level 
# = 1.5 mW-cm~^. For clarity a grid has been superimposed on the figures. Grid size (a) approximately 20 x 20mm. 
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VII. CONCLUSION 

We have identified glider gun-like structures in experimental and numerical studies of the light-sensitive BZ reaction. 

In experimental systems with a heterogeneous network of excitable and non-excitable cells imposed, spirals rotating 
within excitable cells (glider guns) formed. These stable glider guns emit streams of wave fragments into all neighboring 
excitable cells. In numerical studies where the excitability of the network was high we observed similar stable glider 
gun formation. However, in numerical studies we lowered the overall excitability of the network such that certain 
junction transitions become unfavorable and glider guns with fewer stable streams were observed. Glider guns with 
intermittent firing patterns were also observed. This intermittent firing was due to the complex periodic rotation of 
the spiral within the excitable cell. 

We showed in numerical studies how unstable glider guns could be used to implement colhsion-based gates via 
interactions with other traveling wave fragments. Therefore, theoretical ideas concerning universal computation in 
these systems is closer to being realized experimentally. We were able to manipulate glider streams, for example 
annihilate selected streams and switch periodically between two interacting streams. We were also able to show that 
glider guns could be formed or annihilated via specific interactions with glider streams from a second gun. We also 
showed examples where glider guns could be used to implement simple memory analogs. These observations highlight 
the potential for computing with such systems. 

To truly embrace the architectureless philosophy of collision-based computing it will be necessary to identify glider 
gun-like structures in homogeneous BZ systems. In experiments using a weakly-excitable BZ system we identified 
unstable spiral waves that spontaneously break up to form continuous streams of fragments. Both this type of glider 
gun and the types previously identified in heterogeneous BZ systems are based on rotating spiral waves. These spiral 
waves periodically lose stability generating streams of excitation wave fragments. This mechanism has strong parallels 
with CA models based on chemical systems where gliders are formed from dynamical instabilities in rotating spiral 
waves. 

Future work will involve a continued search for unstable glider guns in heterogeneous experimental systems. A 
reproducible and controllable method for forming glider gun analogs in homogeneous BZ systems will also be pur- 
sued. This will achieve the goal of having a truly architectureless computational medium with the well documented 
advantages that this provides. In the future, we aim to construct prototypes of arithmetical chips in the BZ medium. 
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